Theory of vortices in hybridized balhstic/diffusive-band superconductors 
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We study the electronic structure in the vicinity of a vortex in a two-band superconductor in 
which the quasiparticle motion is ballistic in one band and diffusive in the other. This study is 
based on a model appropriate for such a case, that we have introduced recently [Tanaka et al., Phys. 
Rev. B 73, 220501(R) (2006)]. We argue that in the two-band superconductor MgB2, such a case 
is realised. Motivated by the experimental findings on MgB2, we assume that superconductivity 
in the diffusive band is "weak," i.e., mostly induced. We examine intriguing features of the order 
parameter, the current density, and the vortex core spectrum in the "strong" ballistic band under 
the influence of hybridization with the "weak" diffusive band. Although the order parameter in the 
diffusive band is induced, the characteristic length scales in the two bands differ due to Coulomb 
interactions. The current density in the vortex core is dominated by the contribution from the 
ballistic band, while outside the core the contribution from the diffusive band can be substantial, 
or even dominating. The current density in the diffusive band has strong temperature dependence, 
exhibiting the Kramer-Pesch effect when hybridization is strong. A particularly interesting feature 
of our model is the possibility of additional bound states near the gap edge in the ballistic band, 
that are prominent in the vortex centre spectra. This contrasts with the single band case, where 
there is no gap-edge bound state in the vortex centre. We find the above-mentioned unique features 
for parameter values relevant for MgB2. 

PACS numbers: 74.20.-z, 74.50.+r, 74.70.Ad, 74.81.-g 



I. INTRODUCTION 

One can learn a great deal about the pairing mech- 
anism of a superconductor from the electronic proper- 
ties in the presence of inhomogeneity, such as vortices 
and impurities. Although multiple-band superconductiv- 
ity was first studied almost 50 years ago,^ the structure 
of a vortex in multiband superconductors has not been 
well understood, especially when impurities are present. 
The best material that has been discovered so far for 
studying multiband superconductivity is MgB2/^ The 
consensus is that superconductivity is driven by electron- 
phonon interactions, and that it can be well described 
by a two-band model, with the "strong" a band (energy 
gap A,^ « 7.2 meV) and the "weak" tt band (A.^^ ~ 2.3 
meV)<^i^i^ii^'^'i^iii'i^'i^ The two energy gaps vanish at a 
common transition temperature Tc (Refs. HflslflGl ) and 
there is evidence of induced superconductivity in the tt 

The one unusual aspect of this material as a multiband 
superconductor is the effect of impurities. The standard 
theory2i2£iiSl tells us that interband scattering by non- 
magnetic impurities should reduce Tc and the gap ratio. 
Mazin and co-workers, ^^-^^ however, have shown theo- 
retically that this docs not apply to MgB2, due to dif- 
ferent symmetries of the a and tt orbitals, and hence 
negligible interband scattering. Moreover, impurities 
or defects, in particular those at Mg-sites which tend 
to occur more easily than at B-sites, affect only the tt 



band strongly. Indeed, many experiments have shown 
that the a and tt bands are essentially in the ballistic 
and diffusive limit, respectively^^'^i^'Se^'^s^^^ 
Even in samples in which the a band is influenced 
by impurities substantially, two-gap superconductivity 
is retained^i^^'^^i^^^^^iiadi The question of 
whether or not the two gaps merge in dirty samples has 
not yet been settled^^iiiii^ii^ 

Due to induced superconductivity and possibly also 
owing to these peculiar effects of impurities, the elec- 
tronic structure around a vortex in MgB2 has been found 
to exhibit intriguing properties. Eskildsen et al.^ by 
tunneling along the c axis, have probed the vortex core 
structure in the vr band. The local density of states 
(LDOS) was found to be completely flat as a function 
of energy at the vortex core, showing no sign of bound 
states. Also the 'core size' (^50nm) as measured by a de- 
cay length of the zero-bias LDOS turned out to be much 
longer than expected from Hc2 (^lOnm)j ^^i"^'* Moreover, 
the existence of two effective coherence lengths in MgB2 
has been suggested by the ^SR measurement of a vor- 
tex lattice4^ There has been, however, no experiment 
directly probing the a band in the vortex state, and the 
electronic structure around a vortex in the u band is yet 
to be determined. 

Theoretically, the vortex structure in a two-band su- 
perconductor has been studied in terms of two clean 
bands^^ and two dirty bands. Neither of these mod- 
els, however, applies to many MgB2 samples, which have 
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the clean a and dirty tt bands. Recently, we have for- 
mulated a unique model for a multiband superconductor 
with both a clean and a dirty band.— We have studied 
the effects of induced superconductivity and impurities 
in the weak diffusive band on the electronic properties in 
the strong ballistic band around a vortex. A particularly 
intriguing feature found in this model is the possibility 
of bound states near the gap edge in the ballistic band in 
the vortex core, in addition to the well-known Caroli-de 
Gennes-Matricon bound states.^ Such bound states do 
not exist in a single ballistic band, and they arise solely 
from coupling to the diffusive band. Our model has also 
been applied to study the Kramer-Pesch effecli^ in cou- 
pled clean and dirty bandsi^ It has been found that the 
Kramer-Pesch effect is induced in the dirty band, which 
is absent when there is no coupling with the clean band. 
Thus, hybridization of ballistic and diffusive bands can 
lead to unusual properties of a multiband superconduc- 
tor, and the vortex core structure is an example in which 
the effects of induced superconductivity and impurities 
manifest clearly. 

In this work, we make an extensive study of the 
electronic structure in the vicinity of a vortex in a strong 
ballistic band, under the influence of hybridization with 
a weak diffusive band. Our model is based on coupled 
Eilenberger and Usadel equations, which are solved 
directly and numerically. We assume that interband 
scattering by impurities is negligible and that the two 
bands are coupled only by the pairing interaction, as 
justified for typical MgB2 samples. Unique features 
of the order parameter, the current density, and the 
vortex core spectrum are examined in detail. In par- 
ticular, we study the development of the gap-edge 
bound states as various physical parameters are varied. 
Although this work has been motivated by the impurity 

^g?^,^,^ jl6. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. 34. 35. 36. 37.38. 39.40. 41 

and the STM measurement of the vortex state in 
MgB2,^^'^'* our model is general and applicable to any 
superconductor in which ballistic and diffusive bands 
are coupled mainly by the pairing interaction. 

The paper is organized as follows. The formulation 
and computational details are given in Sees. |TT] and IIIIl 
respectively. Results are presented in Sec. IIVI and they 
are summarised and discussed in Sec. [V] Throughout the 
paper we use units with h = c = I. 



II. THEORETICAL DESCRIPTION 



matrix structure of the propagator in the Nambu-Gor'kov 
particle-hole space. In the ballistic case, the equation of 
motion for g is the Eilenberger equation, and in the diffu- 
sive case the Usadel equation. Our model is appropriate 
for any two-band superconductor with a clean and a dirty 
band (generalisation to several bands is straightforward). 
However, having in mind MgB2, for definiteness we call 
the two bands a and tt bands, respectively. 

In the clean a band, P_Fo-, R) satisfies the Eilen- 
berger equation,—!^ 



era - A^, g„ + ivpa ■ ^Qa = 0, 



(1) 



where \'F(j is the Fermi velocity and iS.„ the (spatially 
varying) order parameter. The three Pauli matrices in 
Nambu-Gor'kov space are denoted by fi, i = 1,2,3, and 
[...,...] denotes the commutator. Throughout this work, 
we ignore the variation of the magnetic field in the vortex 
core, assuming a strongly type-II superconductor (this is 
justified for example for MgB2). 

For the tt band we assume that it is in the diffusive 
limit. In the presence of strong impurity scattering, 
the momentum dependence of the quasiclassical Green 
function is averaged out, and the equation of motion for 
the resulting propagator g^(e,R) reduces to the Usadel 
equation's 

efg - A^, + V-(.g^Vg^) = 6, (2) 

L J TT 

with the diffusion constant tensor D. Both ballistic and 
diffusive propagators are normalized according to^ 



~2 -2 2i 



(3) 



A two-band superconductor with a ballistic and a dif- 
fusive band can exist only if interband scattering by im- 
purities is weak. We neglect in the following interband 
scattering by impurities, and assume that the quasipar- 
ticles in different bands are coupled only through the 
pairing interaction. Self-consistency is achieved through 
the coupled gap equations for the spatially varying order 
parameters in each band. 



A«(R) =^ya0iVF;3.f;3(R), 

/3 



(4) 



where a, /? G {cr, tt}, and Aq = fi Re Aq,— f2 Im A^. The 
coupling matrix Va0 determines the pairing interaction, 
Npp is the Fermi-surface density of states on band /3, and 



Both the ballistic and diffusive limits of superconduc- 
tivity can be described within one unified theory, the qua- 
siclassical theory of superconductivity .^-'^'^^i^^'^'^i^^'^^i^''' 
The central quantity of this theory, containing all the 
physical information, is the quasiclassical Green function, 
or propagator, ^(e, pi?Q., R). Here e is the quasiparticle 
energy measured from the chemical potential, p^Q, the 
quasiparticle momentum on the Fermi surface of band a, 
and R is the spatial coordinate. The hat refers to the 2x2 



^a(R) 



^ (/<T(e,PF^,R))pF^tanh( ^ 



•^.(R) = / |^/-(^'I^)tanh(^ 



(5) 



Here is the upper off-diagonal (1,2) element of the ma- 
trix propagator g^, and is a cutoff energy. The Fermi 
surface average over the cr-band is denoted by (• • • )pp„ . 
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Note that within our model the precise form of the 
Fermi surface in the tt band is not relevant, as in the diffu- 
sive limit all necessary information is contained in the dif- 
fusion constant tensor B. For simplicity we have assumed 
in our calculation an isotropic tensor Uiij — DSij. The 
diffusio n constan t defines the 7r-band coherence length 
= \/D/2ttTc. For the a-band Fermi surface we as- 
sume a cylindrical shape, as motivated by the Fermi sur- 
face of MgB2. This allows us to treat the a band as 
quasi-two-dimensional, Vp^. — vp^Pr, with cylindrical 
coordinates (pr,P<t>^Pz) and the unit vector in direc- 
tion of Pr- We define the coherence length in the a band 
as = vpa II-kTc- It will be used as length unit through- 
out this paper. 

The numerical solution of the (nonlinear) system of 
Eqs. H])-© is greatly simplified by using the Riccati 
parameterization of the Green functions, both for the 
ballistic cas o^^'^^i^°'^^ and for the diffusive casci^ 



ga = 



(6) 



1 - 7a 7a 27q 
1 + 7a7a V 27a 7a 7a " 1 

with 7cr(e,pFCT,R) = 7*(-e*, -p^cr, R) and 7^(e, R) = 
7*(— e*,R). The transport equation for 7cr(e, Pfo-, R-) is 
given by^*59.60^ 

A, + 2e 7, + A;72 + zvj., V7, = 0. (7) 
and for 7Tr(e,R) it is^^ 

A. + 267. + A:7^-*i?fvV-^?^^^^') =0. (8) 

We solve Eqs. (P)-® self-consistently. After self- 
consistency has been achieved for the order parameters, 
the (for the a band, angle-resolved) LDOS in each band 
can be calculated by 

Na{€,ppa,R)/Np„ = - Im g„(e,pFa,R)/7r, 

N^ie,R)/Np^ = - Im5,(e,R)/7r, (9) 

where ga is the upper diagonal (1,1) element of ga- 

The current density around the vortex has contribu- 
tions from both tt and a bands, 

j(R)=j,(R)+j,(R). (10) 

The corresponding expressions are 

rim .g^)pF„ tanh (^^) , 



j.(R) 


/•°° de 


2eNp„ 


i-co 2^ 


j.(R) 




2eNp^ 


TI" J-00 '■ 



where e = — |e| is the electron charge. 

III. COMPUTATIONAL DETAILS 

A. Gap equations 

The first step is to write Eqs. (H]) and ([5]) in a form 
that is independent of the cut-off energy Cc and the cut- 
off dependent interaction matrix Vap ■ We diagonalize the 



interactions by 

Np ^v^p^py^^^ ^ylt^\'^^\ (12) 

where Np = Np„ + Np^^ and — Np^/Np. The band 
index runs over a — (ct, tt) , and thus we need two basis 
functions for a complete system, denoted by = 0, 1. We 
expand the order parameter and the anomalous Green 
functions in terms of the eigenvectors 



y - 



yi'K 



(13) 



The eigenvectors are orthonormal. We also introduce the 
vectors 



A = 

Using the expansion 



71 (J J' fj 



(fe) 



A= ^^''^ y ' ^ 



fe=0,l 



k=0,l 



(fe) 



(14) 



(15) 



with A^^) =3^ A and T''^^ =y T, the gap equations 
are given by 



A('=)(R) = A('^)j^('^)(R) 



(16) 



Without restriction we can choose the largest eigenvalue 
to be A^*'). It can be eliminated together with the cutoff 
frequency in favor of the transition temperature between 
the normal and superconducting states, T^. Although 
both coupling constants A^"-* and A'-^-' are cutoff depen- 
dent, the parameter A defined by 



A 



aWa(i) 

A(o) - A(i) 



(17) 



is cutoff independent. We parameterize the two coupling 
constants by the two cutoff independent quantities Tc 
and A. 

Near Tc, the subdominant order parameter A^^-* is of 
higher order in (T^ — T) than A^^^ is. The bulk gaps on 
the two Fermi surface sheets are thus determined near 
by A*^*^) only, and are given by 



y. 



(0) 



y-; 



(0) 



(18) 



We define the ratio of the bulk gaps on the two bands 
near Tc, 



A„ 




(19) 



We also introduce the notation n := ^Jriajn^. Then the 
matrix of eigenvectors is given by 



y'^' y^^ 



-p 

n 



(20) 
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B. Homogeneous solutions 

In the homogeneous case, at zero temperature 
hm Ta = In — — - , 

T^O \AJ 



(21) 



where Ec is the usual BCS cutoff frequency of the order 
of the Debye frequency. Using the BCS formula 



In 



2£c 



1 



7 



(22) 



(7 « 0.577 is Euler's constant), we obtain the zero- 
temperature gap equations 

7 + h.^) AO .-lV3^W3;«%/^AO (23) 



TlTr 



where A° denotes the bulk order parameter at zero tem- 
perature on band a. From this we can derive an equation 
for the ratio po = A°/A° (assuming that the gap ratio 
p is real). 



A- 



Po 



■ln|po| 



P- Po 

'n? + p2 71" + Pop 
and for the dominant zero-temperature gap 



(24) 



(25) 



Equation (|24p is a quadratic equation for /?, which is read- 
ily solved analytically in terms of poj a-nd A: 



2n2po(l + Aln|po|) 



(n2 - pI) + ^(n2 +p^)2 _ (2npoAln|po|)^ 



(26) 



The four material parameters po: A, and com- 
pletely specify the bulk behaviour of the system. We 
illustrate in Fig. [1] how the dominant bulk gap to ra- 
tio at zero temperature is affected by poi f^Tr/^^cr = 1/"-^, 
and A. In (a), A^jTc determined by Eqs. and ([26]) 
is shown as a function of poi for n-,^ — n„ and various 
values of A. In (b) we compare A^/Tc as a function of 
riTr/nrj for several values of po, and for A = —0.1 and 
0.1. It can be seen in Fig. [IJa) that for a fixed n^jn^^ 
A^^/Tc exhibits nonlinear behaviour as a function of pa 
depending on the value of A, while the BCS value (1.764) 
is recovered at po = and 1. For A < this ratio is al- 
ways larger than the BCS value, and it is always larger 
for A < than for A > 0, for fixed po and n^/rio-- When 
A > 0, A^/Tc can be smaller than the BCS value for a 
certain range of po- The effect becomes even more pro- 
nounced for larger nT^jn^. For a fixed poj as shown in 
Fig. [T]Jb) , A^ /Tc increases almost linearly with increas- 
ing n^/ricr, with larger slope for larger po- Here again, 
the enhancement of Aj^ /T^ is always stronger for A < 
than for A > 0. For nTr/n„ = there is only one band 
and the ratio reduces to the BCS value. 




FIG. 1: (Color online) The bulk \/\„\/Tc in the zero- 
temperature limit as given by Eqs. (|25p and (I26p : (a) as a 
function of po ~ AJJ/Act for — and various values of 
A, and (b) as a function of n^/no- for various sets of po and 
A. For a fixed nn/na, \Aa\/Tc shows nonlinear behaviour as 
a function of po, where the BCS value is recovered at po = 
and 1. The bulk gap is enhanced as n-n /ria increases, and this 
effect is larger for larger po and for smaller A. 



In fact, the enhancement of the bulk zero-temperature 
gap by stronger Coulomb interactions in the subdomi- 
nant channel can be estimated for small |A| analytically 
from Eqs. (pS)) and which give 



dk Tc 



A=0 



A^ 



A=0 



npo In I Pol 
n2 + pI 



(27) 



and A°(A) = poA°(A). More generally, at low T for 
given Po and n, the homogeneous gap over ratio in- 
creases with decreasing A in both bands. The slope at 
Po = 1 in Fig. [IJa) is independent of A and is given by 
— 7re^'''/(l -t- 71,2). An expansion for small po or for small 
nTi/n^ gives for both cases 



A^ 



: 7re 



-7 



'l-^(l + Aln|po|)ln|po|pg' 



(28) 



We now discuss the relation between the gap ratios 
near T^, given by p, and at T = 0, given by po. When 
A = 0, we have p = po. For small |A| or small po/n we 
obtain from Eq. (|26p the approximate relation 



pwpo + Apo ln|po| 



(29) 



For I Pol < 1 it follows that when the Coulomb repulsion 
dominates in the second channel (A < 0), the magnitude 
of the gap ratio at Tc is increased with respect to that 
at zero temperature, |p| > |po|; for effective attraction in 
the subdominant channel (A > 0), |p| < |po|. 

In Fig. [5] we show the solutions of the homogeneous 
gap equations for both a and tt bands as a function of 
temperature for various values of A, for fixed po = 0.3 
and n — 1. In the left-hand and right-hand panels of 
Fig. [2] the order parameters in the two bands and the 
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FIG. 2: (Color online) Energy gaps A^(T) and A^(T) (left- 
hand panel) and the ratio Ati{T) / lS.a{T) (right-hand panel) 
in the homogeneous case as a function of temperature T for 
various values of A, for fixed po = 0.3 and n = yjua juT^ — 1. 
There is a second transition temperature associated with the 
TT band for A > Ac (see text). 



ratio of the two are plotted, respectively, as a function 
of temperature. As discussed above, when A = (thick 
lines), the ratio |A7r(T)/Acr(T)| stays constant and equal 
to p = Pa, while for A < (> 0) p is larger (smaller) than 
Po- 

If the order parameter in the tt band has a zero slope 
at Tc, i.e., p — (see, e.g., the lowest curves in Fig. [5]), 
then for positive A there is a second superconducting 
transition^ associated with the tt band at temperature 



TcAp = 0) = T,e-i/A (A > 0). 



(30) 



In this case the order parameter in the tt band stays zero 
above Tctt', below Tctt both order parameters are nonzero, 
and the zero-temperature gap ratio is po = e~^l ^ so that 
the relation Tc-kIT,, = pQ holds. For negative A there is no 
second transition in the tt band, and superconductivity is 
purely induced for all temperatures. For small absolute 
values of p, the order parameter in the tt band is given 
approximately by 



A.(r) 
A.(r) 



1 - A 



In. 



i\p\ « !)■ 



(31) 



This equation is correct for T > T^tt as long as 
\A.^{T)/Aa{T)\ stays small. For p the order param- 
eter in the tt band is nonzero in the entire temperature 
range [see Fig. [21 a)], being purely induced for A < 0, and 
for T > Tc-;, when A > 0. It follows from Eq. jSl]) that 
for A < 0, |A7r(T)/Ao.(T)| is reduced from \p\ as temper- 
ature decreases, while for A > it increases from \p\ for 
<T < Tc [see also Fig. ^h)]. 
For a given po, the critical A, above which a sec- 
ond transition in the tt band can be observed, is given 
by Ac = —1/ In I Pol- In the case of Fig. [2] this yields 
Ac ~ 0.8306. We find that for A > Ac there is a tt phase 
difference between the order parameters in the two bands 
for T > Tc-K (the gap ratio becomes negative), and for 




FIG. 3: (Color online) The A°/Tc (top) and l^l/T^ (bottom) 
as a function of the gap ratio at Tc, p — limT-.Tc Att/Act, 
obtained from Eqs. and ((26]), for A < [(a) and (c)] 

and A > [(b) and (d)]. While for -1 < A < there is a 
unique solution, for A > 0, for sufficiently small p there are 
three solutions. This indicates possible first-order transitions 
as a function of temperature, if A for a given p is larger than 
a critical value. 



T < TcTT a more complicated picture arises, with three 
solutions that compete, and possible first-order phase 
transitions. This is demonstrated in Fig. [31 in which 
the zero-temperature bulk gap to Tc ratio in the a (top) 
and TT (bottom) band is plotted as a function of the gap 
ratio at Tc, p — lim7-_^7-^ A^/A^-, obtained from the so- 
lution of Eqs. and (|26p. Results are shown for the 
effective subdominant coupling A ranging from repulsive 
(left) to attractive (right). For — 1 < A < there is a 
unique solution for the bulk gaps at zero temperature for 
a given p (left-hand panels). However, as can be seen in 
the right-hand panels, when A > 0, there are three pos- 
sible solutions to the gap equations for sufficiently small 
p. As these zero-temperature solutions are associated 
with the same single solution near Tc, this indicates that 
first-order phase transitions can occur as temperature is 
varied. This behaviour is not the topic of this paper, 
however, and we assume for the remainder of the paper 
that A < Ac for a given po (or pQ > p^ — e^^l^ for a 
given A). 



C. Choice of material parameters 

For MgB2, if estimated from the ah initio values 
of the electron-phonon coupling strengths^iSiiiS and the 
Coulomb pseudopotentials,~'S3- A varies depending on 
which calculations are used, while some calculations yield 
very small A. But using Eq. one can fix A from the 
observed values of po and p - although for the latter, some 
data points have relative ly larg e error bars. From the ex- 
periments of Refs. H, ir^fisllla . and HI, we find roughly 
—0.3 < A < 0.3. Among these, the measurement of Ref.[^ 
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has a data point at the highest temperature, T ~ 0.96Tc, 
and A estimated using this point is close to zero. There- 
fore, we assume that for MgB2, A is small and can be 
either positive or negative. 

The strength of induced superconductivity in the tt 
band is specified in our model by A and po, the latter of 
which is about 0.3 in MgB2. Results are presented mainly 
for pq = 0.1,0.3,0.5. The n-^/na in MgB2 has been de- 
termined experimentally to be n-^/na ~ 1 — 1.2Ji^i^ 
We show results for nj^/n^ — 1 and 2, and discuss the 
effects of the 7r-band density of states. 

Finally, there is another material parameter in our 
model, that is the ratio of the coherence lengths in the 
two bands, £,Tr/£,a-- Note in this regard that in defining 
the TT-band coherence length = •\/D/27r Tc, we have 
used the energy scale Tc rather than A^. This is mo- 
tivated by the fact that, as discussed above, supercon- 
ductivity in the tt band is assumed to be mostly induced 
by the a band. Consequently, for the variation of the tt- 
band order parameter determined by the self-consistency 
equation, Eq. (jH), the energy scale Tc is relevant. How- 
ever, the quasiparticle motion in the tt band is governed 
by the Eilenberger transport equation supplemented with 
an appropriate impurity self-energy E^^^ of the order of 
I/ttt, where is the 7r-band quasiparticle lifetime due to 
elastic intraband impurity scattering. To determine un- 
der which conditions the diffusive approximation can be 
used in the tt band, the relative size of the 7r-band gap, 
Att, and the tt band impurity self-energy, E^^^, matters. 
Thus, the condition for the tt band to be diffusive so that 
the Eilenberger equation supplemented with E^^^ reduces 
to the Usadel equation is 



l/r^ > A^. 



(32) 



In practise it turns out that the diffusive approximation 
applies quite well already when is only larger, not 
much larger, than A^r. Using the above definitions for ^o- 



and and with the diffusion constant D 
yields 



3 ' ^ "Fit 



, this 



3A, 



(33) 



Strong electron-phonon coupling in the a band renor- 
malizes the Fermi velocity as vpajZ^^ where Z is 
the renormalization factor and is about 2 in the zero- 
frequency limiti^ Taking the ah initio Fermi velocities, 
VF-n/vpa = 5.8/4.4,~ as the unrenormalized values, and 
using A^/Tc « 0.6 as found in several experiments, we es- 
timate ^7r/C<T < 5 for MgB2. As discussed in Sec. IIVC H 
an estimate from the experiment of Ref. [13 yields a value 
between 1 and 3 for MgB2. To keep the discussion gen- 
eral, the effects of the coherence length ratio on the vor- 
tex core structure are illustrated for ^tt/^o- = 1, 2, 3, 5. 

Throughout this work, cylindrical symmetry is as- 
sumed for the two-dimensional cr band, and an isotropic 
diffusion tensor for the three-dimensional tt band. Up to 
this point, the energy gaps (homogeneous order param- 
eters) in the two bands have been denoted as Ajr and 



1.5- 



-,'^0.5 




FIG. 4: (Color online) (a) Order parameters |A^(a;)| and 
|Ao-(2:)| and (b) the ratio |A^(2;)/Aa-(a;)| as a function of co- 
ordinate X along the path y = G for various values of A, for 
Cir/C<T = 3, po = 0.3, r/Tc = 0.1, n = 1. The Coulomb re- 
pulsion (which reduces A) renormalizes the recovery lengths 
of the order parameters in the two bands differently. 



Ag-. In the following, A^r and A„ represent the spatially 
varying order parameters. 



IV. RESULTS 
A. Order parameter 

1. Effects of tfie Coulomb interaction in the diffusive band 

Negative A means that the effective Coulomb interac- 
tion dominates over the effective pairing interaction in 
the second pairing channel. In this case, Coulomb in- 
teractions are responsible for the fact that the tt band 
cannot maintain superconductivity on its own, and it 
superconducts only due to its proximity to the super- 
conducting a band. We first examine the effects of the 
Coulomb interaction on the order parameters. We con- 
sider a vortex that extends in z direction, with the vortex 
centre situated at a; = ?/ = for each z. In Fig. 21 we 
show the order parameter magnitudes as a function of 
coordinate x along the path y — Q for three values of 
A, appropriate for attractive, repulsive and no effective 
interaction in the second pairing channel. Figure [H^a) 
shows |A7r(x)| and |A(j(a;)|, and Fig. IDJb) their ratio. The 
Coulomb repulsion renormalizes the recovery lengths of 
the order parameters (the characteristic length scale over 
which the order parameter recovers to the bulk value) in 
both bands. Although superconductivity in the tt band 
is mostly induced by that in the a band, it is clear from 
Fig. d] that the length scales in the two bands can differ 
considerably if A ^ 0. For A = there is only one length 
scale and hence the same recovery length in both bands. 
When A < (> 0), |A^(a;)| is enhanced (suppressed) 
around the vortex core, resulting in a shorter (longer) re- 
covery length compared to the case with A = 0. On the 
other hand, the sign of A has opposite effects on |Ao-(x)|, 
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FIG. 5: (Color online) (a) Order parameters |A^(a;)| and 
|Act(2;)| and (b) the ratio \/S.Tr{x) / Aa(x)\ as a function of co- 
ordinate X along the path y = for various density of states 
ratios nT^jn^, and fixed Cir/Co- = 3, po = 0.3, T/Tc = 0.1, 
A = —0.1. The TT-band density of states renormalizes the 
recovery lengths of the order parameters, and the bulk mag- 
nitudes of the order parameters are larger for larger n-n jric. 



although these effects are relatively small for the parame- 
ter set for Fig. [31 In fact, A^- (its length scale and magni- 
tude) can be affected by the 7r-band Coulomb interaction 
significantly for larger po and ^^/^^-^ 

As demonstrated in Fig. [1] the bulk gap over Tc ratios 
in both bands are enhanced by stronger Coulomb repul- 
sion in the diffusive band, an effect also seen in Fig. UJ^a) . 
As shown in Eq. (j27p , this renormalization of the magni- 
tude of the bulk gap has the same sign for both bands, 
and is very different from the renormalization in the core 
region seen in Fig. Ufa), that has opposite sign for the 
two order parameters. The latter effect can be seen most 
dramatically when plotting the ratio |A^(a;)/Acr(a;)|, as 
we do in Fig. [Hb) . For A 7^ the ratio increases by up 
to a factor of 2 for A ~ —0.1, and decreases by about the 
same factor for A = 0.2. When A = 0, \1\t,{x) j I\^{x)\ 
stays constant and equal to po = P- This results from 
the fact that in this case A*^^-* = in Eq. HI]), leading to 
A(^)(j;) = so that Eq. HH]) holds for aU x. 



2. Effects of the density of states in the diffusive band 

Figure [5] demonstrates the effects of the Fermi-surface 
density of states in the diffusive band on the order pa- 
rameters for a typical set of material parameters and for 
T/Tc ~ 0.1. In Fig. El^a) we show the order parame- 
ter magnitudes and in (b) their ratio. It can be seen 
in Fig. EJa) that the bulk gap over Tc ratios increase 
with increasing density of states ratio nT^jria, a behavior 
consistent with Fig. [TJb) and with Eqs. and (^5)1 . 
From Fig. [SKa) we see that a larger tItt/^^ct also enlarges 
the core area in both bands. At higher temperature 
(not shown), n^/rto- has less effect on the order param- 
eters. The effects of the 7r-band Coulomb interaction 
described in the last section manifest more drastically 
in both bands for larger n-^/na. When comparing the 
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FIG. 6: (Color online) (a) Order parameter in the tt band 
|A.n-(a::)| as a function of coordinate x along the path y = 
for various values of po and ^ir/Co-, for T/Tc = 0.1 and 
rin/n^ = 1; for (a) A = -0.1 and (b) A = 0.1. For A > 0, 
jA-^(x)| is depleted around the vortex core for larger £,-n/(,a, 
and this effect is larger for larger po- When the Coulomb 
repulsion dominates (A < 0), for small po (e.g., 0.1), |A,r(a;)| 
is enhanced near the vortex centre for larger 5,r/Co-. 



shapes of the order parameter profiles in the a band and 
the TT band in Fig. [5l^a) , there is a qualitative difference 
between the cases n^/ng- = 0.5 and n-^/ucr = 2. Never- 
theless, Fig. [DJb) shows that despite the different shapes 
of the order parameter profiles in the two bands for dif- 
ferent values of nj^/n^, the spatial variation of the ratio 
\At^{x)/ A^{x)\ is almost independent of n-^ln„. This is 
for A a nontrivial and rather interesting finding. 



3. Effects of the coherence length and gap ratios 

In Fig. [6] we illustrate how the order parameter profile 
in the tt band changes as po and ^jr/Co- are varied. We 
compare the two cases of (a) A = —0.1 and (b) A = 0.1, 
leaving T/T^ = 0.1 and n^/no- = 1 fixed. As can be 
seen in Fig. [HJb), for A > 0, |A^(a;)| is depleted in the 
core area as ^tt/Cct increases, and this effect is larger for 
larger po- On the other hand, when the Coulomb re- 
pulsion is dominant in the tt band (A < 0), as shown 
in Fig. ini^a), the profile of |A^(a;)| changes in a peculiar 
way as ^t/'Cct varies. For relatively small po [e.g., 0.1 in 
Fig.mja)], |A^(a;)| is enhanced around the vortex core for 
larger ^tt/^ct, while for po = 0.3 |A7r(a;)| hardly changes 
for different values of £,tt/£,(j- For larger po(> 0.3) (not 
shown), however, the core region becomes larger with in- 
creasing ^7r/^CT - in a similar way as for A > 0, though 
not so drastically. This lets us conclude that for the case 
Po ~ 0.3 in Fig. [HUa) a cancellation between two oppo- 
site tendencies is at work, that renders |A^(a::)| seemingly 
insensitive to the parameter ^jr/Co-- to |Ao-(a;)|, simi- 
larly as |A^(a;)| in Fig. [6l[b), the core area is enlarged for 
larger £,T^/(,a, more noticeably for larger po, and for both 
positive and negative A [see Fig. [Htb) below] . 

As discussed above, even though superconductivity is 
induced in the tt band due to coupling with the a band. 
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FIG. 7: (Color online) (a) The ratio |A^(2;)/A^(a;)| as a func- 
tion of coordinate x along the path y = for various values 
of /9o and (^■n/£,a, for T/Tc = 0.1 and n-njua = 1; for (a) 
A = —0.1 and (b) A = 0.1. For negative (positive) A, the 
ratio is larger (smaller) than po around the vortex core, and 
the area of the enhancement (depletion) is larger for larger 
C^r/^CT. As Po increases, the ratio changes more drastically as 
X approaches zero. 



the characteristic length scale of the order parameter can 
be quite different in the two bands. This is demonstrated 
in Fig. [7] for various sets of po and ^Tr/fcr- Here the order 
parameter ratio |A^(x)/Ao.(x)| is plotted as a function 
of X. In agreement with our discussion of Fig. IDJb), this 
ratio is in the vortex core larger than po for negative, 
and smaller than po for positive A. When A < 0, despite 
the peculiar change of the |A7r(a;)| profile as a function 
of ^tt/Co- depending on the value of po, |A^(a;)/Ao-(a;)| is 
always larger than po in the core area, as can be seen 
in Fig. [Tl^a). The main conclusion of Fig. [7] is that the 
enhancement for A < and depletion for A > of 
\1\t^{x) I t^cr{,x)\ is increasing with ^tt/Cct and pp. 

In Fig. [5] we compare the influence of the material pa- 
rameters Po and ^7r/C<T on the order parameter profiles 
in the vortex core. In Fig. [SJa) we vary po for fixed 
= 3, and in (b) we vary ^tt/^ct for fixed po = 0.3. 
As shown in Fig. [U in general the bulk gap to ratio in 
the tr band changes in a nonmonotonous way as a func- 
tion of pq. In the case of the parameter set for Fig. [5fa), 
however, the bulk |Ao-|/Tc increases with increasing po 
for the values of po shown (for a fixed ^tt/Co-)- At the same 
time, as the coupling of the two bands becomes stronger, 
the recovery length of |Ac,| increases [Fig. [SJa)] and the 
length scales in the two bands are more different [see also 
Fig.[7I[a)]. For a fixed po the cr-band order parameter has 
a larger core area as ^tt/'^ct increases [Fig.[51Jb)]. The |Acr| 
profile changes in a similar way with varying ^tt/^o- when 
A > 0. Also seen in Fig. [8{b) is that the insensitivity of 
the [A^l profile due to the cancellation effect for A < 
discussed above, combined with the enhanced depletion 
of I Act I, yields a seemingly paradoxical picture, namely 
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FIG. 8: (Color online) Order parameters jA-^(3::)| and |ACT(a;)| 
as a function of coordinate x along the path j/ = for T/Tc = 
0.1, A = —0.1, nT^jria = 1; for (a) various values of po for 
£,Tv/£,<y = 3 and (b) different values of ^ir/Co- for po = 0.3. For 
a fixed ^tt/Cct, as po increases, the core area is enhanced in 
both bands. For a fixed po, for larger ^tt/Cct the core area is 
enlarged in the cr band. 
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FIG. 9: (Color online) The ratio |A^(a;)/Ao^(a;)| as a function 
of position x along the path i/ = for Ctt/^o- = 3 and n^/wo- = 
1, for various values of po and T /Tc\ for (a) A = —0.1 and 
(b) A = 0.1. The difference between the length scales in the 
two bands and its temperature dependence are enhanced as 
Po increases, especially for A < 0. As temperature is raised, 
the difference between the recovery lengths of the two order 
parameters becomes smaller. 



that an increase in (for fixed results in a larger 
core area in the a band instead of the tt band. 



4-. Temperature dependence 

The ratio \1^t^{x) / /!^a{x)\ as a function of coordinate x 
is plotted in Fig. [9] for various values of po and T/Tc, for 
(a) A = —0.1 and (b) A = 0.1. The enhancement (deple- 
tion) and its temperature dependence of this ratio in the 
core area are larger for larger po for negative (positive) 
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FIG. 10: (Color online) Order parameters (a) |ACT(a;)| and (b) 
|A,r(2;)| as a function of x along the path y = for po = 0.3, 
^tt/Cct ~ 3, n-n/ria = 1, A = —0.1. The 'weak' tt band makes 
the temperature dependence of the lAo-j profile more drastic 
than in the single-band case. 



A. This effect is considerable when the Coulomb repul- 
sion dominates in the tt band [Fig.[5I^a)]. As temperature 
becomes higher, however, the difference between the re- 
covery lengths of the two order parameters is reduced, as 
can be seen in Fig. [5] (see also Fig. [TU]). 

Figure [TU] shows the temperature dependence of (a) 
I Act (x) I and (b) |A7r(a;)| for a typical set of material pa- 
rameters. By coupling to the weak tt band, the sup- 
pression of I Act I and the widening of its core region with 
increasing temperature are more drastic compared with 
the single-band case. The 7r-band order parameter is less 
affected by temperature for low enough T/Tc- For higher 
temperature, however, the profile of the order parameters 
in the two bands are similar. 

The temperature dependence of the core area of the 
order parameter can be characterised by the vortex 'core 
size' defined by^ 



dA{r = 0) 



dr 



A(r = 00) 



(34) 



where r is the radial coordinate measured from the vortex 
centre. In a clean single-band superconductor, the order 
parameter exhibits the Kramer-Pesch (KP) effect,— i.e., 
shrinkage of the vortex core as T is lowered, approaching 
zero in the zeo-temperature Iimit^i^i22ii21. In an s-wave 
superconductor, however, the core size as defined above 
saturates as temperature approaches zero, when there 
are nonmagnetic impurities. ^"'^^ This can be understood 
by the broadening of the vortex core bound states that 
carry the supercurrent around the vortex centre. This 
broadening removes the singular behaviour in the spatial 
variation of the order parameter in the vortex centre, 
and the vortex core shrinking ceases when ksT becomes 
smaller than the energy width of the zero-energy bound 
states situated at the vortex centre (in quasiclassical ap- 
proximation one can neglect the small splitting of these 
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FIG. 11: (Color online) The core size defined by Eq. §^ 
in the two bands for po — 0.3, Ct/Co- = 3, and n-nln^ = 
1, for (a) A = 0.1 and (b) A = -0.1. The Kramer-Pesch 
effect is induced in the tt band, although absent in a single 
diffusive band (dashed-dotted curves). As seen from (a), the 
core radius is larger in the tt band than that in the a band 
when A > 0. However, the opposite is true for A < 0, as 
shown in (b). 



bound states due to the Caroli-deGennes minigap). 

Recently, our model of coupled clean and dirty bandsi 
has been applied to study the KP effect in a two-band 
superconductor, and it has been shown that, by coupling 
to the clean band, the KP effect is induced in the dirty 
band.^° The calculation in Ref. [sol corresponds to positive 
A in our formulation. 

In Fig. [TT] we plot (,cl(,a in the a and tt bands as 
a function of temperature for po — 0.3, ^jr/^o- = 3, 
and nT,/n^ = 1, for (a) A = 0.1 and (b) A = -0.1. 
Points are data obtained by self-consistent calculation, 
and curves are guides to the eye. Single-band results are 
also shown for clean (dotted curves) and dirty (dashed- 
dotted curves) superconductors, which exhibit the core 
shrinkage and saturation as T — > 0, respectively. In the 
two-band case, the core size in the dirty tt band (dashed 
curves) tends to zero in the zero-temperature limit, as 
the core size in the clean a band (full curves) does. For 
A > 0, shown in Fig. [TTT a). the core size in the tt band 
is always larger than the core size in the a band. There 
is the KP effect in the tt band also when the Coulomb 
repulsion is dominant [see Fig. [TT] (b)] . However, here the 
core size defined by Eq. is always smaller in the tt 
band than in the a band. Surprisingly, with dominating 
Coulomb interactions, the KP effect is better developed 
rather in the diffusive band than in the ballistic band. 
Only at very low temperature the core size extrapolates 
to zero in the a band, whereas this happens in the tt 
band already below T ~ 0.5Tc. We suggest that the KP 
effect can be used to study the strength of Coulomb inter- 
actions in the effective pair interaction matrix. Finally, 
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FIG. 12: (Color online) The 7r-band LDOS iV,r(e,r) as a func- 
tion of energy e at various distances r from the vortex centre 
(r/^o- from to 40 with an increment of 5). The parameter 
values are ^tt/Co- = 1, Po = 0.3, T/Tc = 0.1, n-,tlna = 1, and 
A = 0. The 7V^(e, r) shows no sign of localized states in the 
vortex core, consistent with the experiment on MgB2 (Ref.— ) . 
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FIG. 13: (Color online) The zero-bias LDOS N^{e = 0, r) as a 
function of distance r from the vortex centre for n-,^ jua = 1; 
(a) pa = 0.3, T/Tc = 0.1, A = —0.1 for various values of 
?,r/C<T; (b) i-^lia = 3, T/Tc = 0.1, A = -0.1 for various val- 
ues of po; (c) po = 0.3, Ctt/^ct — 3, T/Tc = 0.1 for various 
values of A; (d) po = 0.3, €,-n/(,(j = 3, A = —0.1 for vari- 
ous t emp eratures. By comparing with the experimental data 
(Ref. [rn) we estimate £,^/£,cr ~ 2 for MgB2, taking po — 0.3. 



when A = (not shown), the order parameters in the 
two bands have the same length scale and the same ^c- 



B. Spectral properties of diffusive band 

In Fig. [12] we show an example for the LDOS iV,r(e, r) 
as obtained from Eq. ([9]) in the diffusive tt band as a 
function of energy e at various distances r from the vortex 
centre. At the vortex centre, Nj^^e^r = 0) as a function 
of e is completely flat, showing no sign of localized states. 
This is consistent with the experiment on MgB2,^^ which 
probed the vr-band LDOS by tunnehng along the c axis. 
Far away from the vortex core, the bulk BCS density of 
states is recovered. For a fixed poi the BCS density of 
states is recovered at larger distances for larger Ctt/Co-; 
while for a fixed ^tt/^o-, the LDOS recovers to the bulk 
density of states further away for smaller po- 

Experimentally the vortex 'core size' obtained from 
tunneling spectroscopy can be defined as a measure of 
the decay of the zero-bias LDOS as one moves away from 
the vortex centre. Note that this definition of the vortex 
core size is very different from that defined in Eq. (I34p 
and shown in Fig. 1111 the two quantities in fact differ by 
a large amount as we discuss in the following. Figure fT3l 
shows the zero-bias LDOS in the tt band NTr{e = 0, r) 
as a function of r for (a) varying ^^/fcr, (b) varying po, 
(c) varying A, and (d) varying T/Tc, for the combina- 
tion of the remaining parameters riT^/n^ = 1, ^tt/Co- = 3, 
Po = 0.3, T/Tc = 0.1, A = -0.1. For a fixed po, the core 
size defined as a decay length of the zero-bias LDOS in- 
creases substantially with increasing ^tt/^ct, as shown in 



Fig.fT^al. On the other hand, for a fixed ■^tt/Co-i the core 
size becomes considerably larger for weaker coupling to 
the a band, i.e., smaller po, as presented in Fig. llSf b). 
It can be seen in Fig. [T3)Jc) that the LDOS profile in the 
TT band is barely changed for different values of A. The 
Fermi-surface density of states n^/no- also has little ef- 
fect on the TT-band LDOS (not shown). As illustrated in 
Fig. [Old), the LDOS is hardly affected by temperature 
over a rather large temperature range; however, the de- 
cay length increases notably for temperatures close to Tc- 
Fig. [T3]is a central result of our calculations, as it can be 
used directly as reference for comparison with scanning 
tunneling experiments, in order to determine either of 
the material parameters ^tt/^o- or po- In particular, the 
parameter ^tt/Co- is otherwise difficult to access experi- 
mentally. 

The decay length of the LDOS can be very different in 
the two bands. The length scale of the zero-bias LDOS 
in the a band is ^o- and thus shorter than that in the tt 
band for small enough po or large enough ^tt/Cct- The 
existence of two coherence length scales in MgB2 has 
been suggested jiii^'^ i.e., a 'core size' (^50 nm) much 
larger than that expected from Hc2 (^10 nm). Theoreti- 
cally, two apparent length scales in the LDOS have been 
found also in the case of two clean bands'*^ and two dirty 
bands4ii2^ The difference in the decay lengths comes 
from the fact that, although superconductivity is induced 
in the tt band, quasiparticlc motion is governed by the 
effective co herence length resulting from the Us adel equa - 
tion dH), ^/D/2A^ (that is longer than = .JD/2ttTc). 

Our model of a clean and a dirty band is suitable for 
describing MgB2, and the values of poj n-j^/na and A for 
Fig- ESIa) are appropriate for modelling MgB2. Using 
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FIG. 14: (Color online) The cr-band LDOS N„{e, r) as a func- 
tion of energy e at various distances r from the vortex centre 
('"/Co- from to 10 with an increment of 1). The parameter 
values are Ct/C"- = 3, po = 0.3, T/Tc — 0.1, rin/ncr = 1, 
and A = —0.1. Coupling to a diffusive band results in bound 
states at the gap edge in the clean band, in addition to the 
well-known Caroli-de Gennes-Matricon bound states. 



the ab initio value^- of vpa- and including the renormal- 
ization factor as mentioned above, £,a — 6.8nm for MgB2 
(Tc = 39 K). By comparing Fig. fTST a) with the zero-bias 
LDOS measured in the experiment of Ref fl7i [Fig. 3(a) in 
Ref. [13], for this parameter set, we find Cn/^a is about 2 
for MgB2. As discussed in detail in the next section, there 
can be bound states near the gap edge in the 'strong' 
ballistic a band, which arise from coupling to the 'weak' 
diffusive tt band. Such bound states exist for parameter 
values relevant for MgB2, e.g., for ^tt/Co- — Po — 0.3, 
and for various temperatures, n„/n„, and A [see, e.g.. 
Fig. [IS]. 



C. Spectral properties of ballistic band 

1. Total LDOS 

We turn now to the spectral properties of the vor- 
tex core originating from the clean a band. We plot 
in Fig. [H the total (angle-averaged) LDOS N„(e,r) = 
{Na[t,TpPcr,r))-pp^ in the ballistic u band as a function 
of energy e at various distances r from the vortex cen- 
tre. The spectra show the well-known Caroli-de Gennes- 
Matricon (CdGM) bound-state bands^i^^ at low ener- 
gies. In the single-band case [see Fig. fTSla)]. at the 
gap edge the spectrum is suppressed and there is nei- 
ther coherence peak nor bound state in the vortex centre 
spectrum,^i25,76 r^^^ 

new feature found in our model is 
the additional bound states at the gap edge, as can be 
seen clearly in Fig. 1141 which arise from coupling to a 
diffusive band. While the CdGM bound states disperse 
strongly as a function of position, the extra bound states 
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FIG. 15: (Color online) The Ncr{e,r) as a function of energy e 
at various distances r from the vortex centre (r/^cr from to 
10 with an increment of 1) for T/Tc = 0.1; (a) the single-band 
case, (b) Ctt/C^ = 5, po = 0.2, (c) ^n/^a = 2, po = 0.3, and 
(d) Ct/Cct = 5, Po — 0.3. In the two-band cases, n^/ria = 1 
and A — —0.1. The gap-edge bound states, which are absent 
in the single-band case, are enhanced (and the number of 
branches increases) for larger ^vr/Co- and po. 



stay near the gap edge as r is varied. The self-consistency 
of the order parameters is essential for studying the pres- 
ence or absence of these gap-edge bound states. 

In Fig.[T3]the total LDOS near the gap edge at various 
distances r is shown at temperature T/Tf, = 0.1, for (a) 
the single-band case, (b) ^^/^.a — 5, po — 0.2, (c) £,Tr/^a = 
2, Po — 0.3, and (d) Ctt/Cct — 5, Po — 0-3. It can be seen 
in Fig. [TSTa) that in the single-band case, there is no 
localized state near the gap edge at r = 0, while a small 
peak develops at the gap edge at large r. This is a bound 
state arising solely from the phase winding around the 
vortex. In Fig. [TSTb'). additional bound states exist at 
the gap edge both in the core region and far away from 
the vortex centre. For a fixed ^Tr/^a, with increasing 
Po, gap-edge bound states are enhanced and more bound 
states can be created [Figs. ITHb) and [TSTd)]. In the 
case of Fig. fWd). there are two bound states close to 
the gap edge in the r = spectrum. For nonzero r. 
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FIG. 16: (Color online) The total LDOS in the a-band at the 
vortex centre, Na{€,r = 0), for energy e close to the energy 
gap, for ^-n/^tT = 5 and T/Tc = 0.1 for different values of A; 
for po = 0.2 (a,b) and 0.5 (c,d) and for n-njua = 1 (a,c) and 
2 (b,d). Bound states exist near the gap edge for A smaller 
than a certain value. 



all bound states broaden into bands due to dispersion of 
the bound states with momentum direction, exhibiting 
the typical l/v^-like behaviour at the band edges. This 
can be seen clearly in Fig. [T3] for the main bound state 
crossing zero energy at r = 0. The additional gap-edge 
bound states also show this behaviour, as can be seen, 
e.g., in Fig. fTSTd) for the main gap-edge bound state at 
small r. For a fixed (more) bound states can exist 
for larger ^,r/Co- [Figs. fTsTc) and [TSFd)]. In the case of 
Fig. fTSfc). gap-edge bound states exist at large r, but 
not in the core region. However, in this case a coherence 
peak can be observed at the gap edge at small r. 



2. Gap-edge bound states 

We now study the detailed development of bound 
states near the gap edge in the a band in terms of the 
spectrum at the vortex centre. We start with the sen- 
sitivity of the gap-edge bound state spectrum to A and 



FIG. 17: (Color online) Same as Fig. \T6\ but for fixed A = 
0, and for (a) various values of ^ir/^o- for po = 0.5 and (b) 
different values of po for ^tt/^ct ~ 5. Gap-edge bound states 
develop as po or ^tt/Co- is increased. 



to nT^/ric. In Fig. [11] we show the cr-band total LDOS 
Ncr{e, r = 0) for energy e close to the energy gap at tem- 
perature T = O.lTc. To emphasize the changes we show 
here examples for a relatively large coherence length ratio 
^^/^^ = 5. The panels (a)-(d) correspond to various fixed 
parameter combinations, and the curves in each panel are 
for different values of A. We show results for po — 0.2 
(a,b) and 0.5 (c,d), and for and for 71^^/11^ — 1 (a,c) and 2 
(b,d). The Coulomb interaction in the tt band also affects 
the gap-edge bound states. It can be seen in Fig. [161 a) 
that, for the given set of parameters, gap-edge bound 
states exist for A = —0.1, but not for A = 0.2. An in- 
crease in n^/no- results in an increase in the energy gap, 
and (more) gap-edge bound states due to widening of the 
core area [Fig. [THT b)]. Increasing po has similar effects. 
It can be seen in Figs. [TCTc.d) that the lowest branch 
shifts considerably to lower energy, while the energy gap 
becomes larger, as A is reduced. 

Next we discuss the dependence on the coherence 
length ratio ^tt/Co- ^-nd on the zero temperature gap ra- 
tio Po- There is a threshold for both parameters, above 
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FIG. 18: (Color online) Same as Fig. [16] but for n^/ua- = 2 
and A = 0, for four sets of other parameters (see text) and 
various temperatures. The energy gap is reduced significantly 
as temperature increases and the energies of the gap-edge 
bound states are shifted together with the energy gap. 



which gap-edge bound states appear. We illustrate this 
in Fig. [17] for T/Tc = 0.1 and A = 0. In Figs. \T7\a.h) 
(.tt/^ct is changed for fixed po = 0.5. Results are for (a) 
n-Tr/na- — 1 and 2 (b). For this gap ratio, gap-edge bound 
states exist for ^tt/Co- > 1 for both ratios of the density of 
states shown. As ^tt/Cct increases, the energy gap changes 
only slightly. On the other hand, the lowest bound-state 
energies are reduced substantially and more branches ap- 
pear. This effect is enhanced for larger riT^/na. For 
Ctt/^ct = 5 and po = 0.5, there are three branches for 
n-^jric = 2. Figures [TWc.d) show the LDOS for various 
values of po for fixed Ctt/Co- = 5 for (c) n^/no- — 1 and 
2 (d). For this parameter set, for po — 0.1 there is no 
gap-edge bound state for riT^jna < 2. As po increases, 
the energy gap is enhanced and more branches appear at 
the gap edge, and this effect is more significant for larger 

Increasing temperature T also enhances the features of 
gap-edge bound states. Figure [TH] shows the LDOS for 
nT^/n„ = 2 and A = for various T. In (a) and (c) we 
compare for ^tt/^o- = 3 the cases pa = 0.3 and po = 0.5, 
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FIG. 19: (Color online) The angle-resolved LDOS Na{e,py = 
0, y) as a function of energy e, for quasiparticles moving in 
the X direction at various positions along the y axis (from 
y — — 18^CT to y = IS^CT with an increment of S^o-)- The 
parameter values are ^vr/^o- ~ 5, T/Tc — 0.5, n-nlric, — 1, 
A = —0.1, with (a) Po = 0.3 and (b) po = 0.5. An increase in 
po results in more gap-edge bound states and enhancement of 
their dispersion. 



and in (b) and (d) we fix po = 0.5 and compare the cases 
^jr/^(T = 1 and £,Tr/^a = 5. As T increases, bound states at 
the gap edge shift lower together with the energy gap and 
increase in strength. The change from T/Tc = 0.1 to 0.5 
in Fig.[TH]is substantial for po — 0.5 and ^■jr/S.a = 5. The 
results in Fig. [TB1[TB1 illustrate that the gap edge bound 
states are a robust feature of our model. 



3. Angle-resolved LDOS 

The bound-state spectrum can be discussed most 
clearly in terms of the angle-resolved LDOS spectra, ob- 
tained from Eq. ([9]). In Fig. [19] the angle-resolved LDOS 
N^{e,Py = 0,y) for quasiparticles moving with impact 
parameter y in the Px direction is shown as a function of 
e at various positions y {x = 0). The two panels are for 
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I ll =5, T/T =0.5, 11 In =1, A=-0.1 




FIG. 20: (Color online) The energy of the bound states as a 
function of y obtained from Fig. 1191 There are two branches 
of gap-edge bound states for po — 0.3. For po = 0.5 there 
is a third branch in the vortex core: the lowest branch is 
extended over large distances and its dispersion is enhanced 
considerably compared with the case for po = 0.3. 



(a) Po = 0.3 and (b) po = 0.5. The position of the bound 
states as a function of impact parameter y is plotted in 
Fig. [20l The CdGM bound-state branch disperses with 
angular momentum, crossing the chemical potential in 
the vortex centre. In Fig. [TOT a) additional bound states 
near the gap edge can be seen, which have only weak 
dispersion. A close inspection reveals that there are two 
branches; for e > 0, the lower branch for roughly y > —6 
and the higher one for y > 0, and similarly for negative 
y for e < 0. The higher branch hardly disperses with y 
[see Fig. [20la)]. As coupling with the tt band becomes 
stronger, the energy gap increases. At the same time, 
the energies of the bound states are lowered (for rela- 
tively small angular momenta), and more bound states 
appear near the gap edge, due to enlargement of the core 
area [Fig. [lUb)] [see also Fig.[20tb)]. For po = 0.5, there 
is a third bound-state branch at the gap edge for small 
|t/|, and the lowest branch extends over large distances, 
with enhanced dispersion compared to that for po = 0.3. 
For a fixed po, increasing ^tt/Co- results in a larger core 
area and has similar effects as increasing po. 



D. Current density 

For a clean single-band superconductor, supercurrents 
around a vortex are strongly coupled to the Andreev 
spectrum discussed in the preceding section, and are 



Po = 0.1 , n^rig = 1 , A = pg = 0.5, n^n„ = 1 , A = 
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FIG. 21: (Color online) Magnitude of the supercurrent den- 
sity in the two bands, j,r and ja, as a function of distance r 
from the vortex centre, for n-nlua — 1 and A = 0, for various 
temperatures. In the left-hand panels, we have po = 0.1, for 
(a) Ct/Cct = 1 ^nd (b) 3. The right-hand panels have po = 0.5, 
for (c) ^tt/^ct ~ 1 and (d) 3. The 7r-band contribution can be 
substantial, or even dominating, and have strong tempera- 
ture dependence far outside the vortex core for large enough 
Po and Cvr/C<T- 



inside the core area carried mainly by these bound 
statesi^iiS In the following we present results for the 
current density calculated with Eqs. pip for our two- 
band model. In Fig. [^T] we show contributions from the 
two bands to the current density around the vortex sep- 
arately, for parameters nj^/na- = 1 and A = 0, and for 
various temperatures. In the left-hand panels, we have 
Po — 0.1, and (a) Ctt/Cct = 1 and (b) 3. For this weak 
coupling between the two bands, the current contribu- 
tion from the tt band is negligible for £,Tr/£,a = 1- For 
Ctt/Cct = 3, though still small, at low temperature (i.e., 
T/Tc = 0.1) the TT-band contribution in the bulk (out- 
side the vortex core, but still well within the penetration 
depth distance from the core) is almost the same as that 
in the a band. As T —^ 0, the peak of jo- approaches 
the vortex centre - this is the KP effect manifest in the 
supercurrent. The current density arising from the in- 
duced superconductivity in the n band is also enhanced 
by decreasing temperature, but the maximum does not 
approach the vortex centre [see Fig. l^lTb)]. In Ref. [l|, 
the same trend of the 7r-band current density as a func- 
tion of temperature was presented for parameter values 
appropriate for MgB2. 

As Po increases, the 7r-band contribution to the current 
density becomes considerable. In the right-hand panels 
of Fig. [51] the current density contributions are plotted 
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for po = 0.5, for (c) Ctt/Cct = 1 and (d) 3. Interestingly, 
as the coupling between the two bands is increased, the 
TT band starts exhibiting the KP effect in the current 
density. The maximum shifts towards the vortex centre 
as temperature is decreased, although the 'core size' as 
defined by the position of the current maximum remains 
larger than that in the a band [see Figs. [2lT c) and (d)]. 

As can be seen in Fig. BTT d) , for relatively large po and 
^^/^^(> 1), the TT-band contribution to the current den- 
sity becomes substantial. Especially away from the vor- 
tex core, where the contribution of the a band is reduced, 
the TT-band contribution can be dominant. Another im- 
portant observation is that the current density in the tt 
band has strong temperature dependence also far outside 
the core. On the contrary, the cr-band current density is 
temperature-dependent only in the core area. Far away 
from the vortex core, we can use in Eq. (|lip the analytical 
solutions to the Eilenberger and Usadel equations for a 
quasihomogeneous system, with constant magnitude and 
constant phase gradient of the order parameters. We 
thus obtain the current-density magnitudes for a large 
distance r from the vortex centre (but small compared 
to the London penetration depth) as 

2„(r) « eNpaVpa^' (35) 

j^{r) « eNF^TTD\A^\- . (36) 
r 

Clearly, the temperature dependence of is dominated 
by that of |A^|. This bulk behaviour of the current den- 
sity in the two bands also explains the dominance of the 
TT-band contribution in the bulk for relatively large pQ and 
^tt/S,(t- For large r, the current-density ratio approaches 

lim _ NF.2TrD\A^\ _ |A^| f C.V 

Thus the TT-band contribution to the current density dom- 
inates when (^^/,^o.)2 > (no-/n^)(|A7r|/Tc)^^ This is rel- 
evant for MgB2, for which this condition roughly reads 
as C^/CtT > 1- 

V. CONCLUSIONS 

In conclusion, we have studied a model recently intro- 
duced by usi for describing a multiband superconductor 
with a ballistic and a diffusive band in terms of coupled 
Eilenberger and Usadel equations. Both equations were 
solved directly and numerically until self-consistency of 
the order parameters was achieved. We have studied the 
effects of induced superconductivity and impurities in the 
weak diffusive (tt) band on the order parameter, the cur- 
rent density, and the spectral properties of the strong 
ballistic (cr) band. A unique feature found in our model 
is the existence of additional bound states at the gap edge 
in the ballistic band, which are absent when there is no 
coupling with the diffusive band. 



Although the two bands are coupled by the pairing 
interaction, the order parameters in the two bands can 
have very different length scales. When the Coulomb 
interaction dominates in the tt band, the order parameter 
in the a band has a longer recovery length than that in 
the TT band. In this case, the tt band exhibits the Kramer- 
Pesch effect with a 'core size' [defined by Eq. ([M]) ] smaller 
than that in the a band. Furthermore, when coupled to 
the diffusive band, the order parameter in the ballistic a 
band is suppressed by temperature more strongly than 
in the single-band case. 

The zero-bias LDOS in the two bands can have very 
different decay lengths, the one in the tt band being larger 
than that in the a band. As the induced superconductiv- 
ity in the tt band becomes stronger, the core area widens 
in the a band and (more) gap-edge bound states ap- 
pear. Moreover, an increase in the Fermi-surface den- 
sity of states in the tt band results in an increase in the 
energy gap and the number of gap-edge bound states in 
the a band. Increasing the Coulomb repulsion in the tt 
band and temperature also enlarges the core area in the 
a band, which results in additional bound states at the 
gap edge. The gap-edge bound states have only weak dis- 
persion. It is thus expected that these bound states are 
affected only weakly by impurity scattering within the a 
band.^^ Results incorporating impurities in the strong a 
band confirm this statement and will be presented in a 
future paper i22 

The supercurrent density is dominated in the vortex 
core by the a-band contribution, and outside the core the 
TT-band contribution can be substantial, or even dominat- 
ing over the cr-band current density. The current density 
in the a band shows the Kramer-Pesch effect in the core 
area. On the other hand, the tt band current density 
has strong temperature dependence also outside the core. 
It exhibits the Kramer-Pesch effect when the coupling 
between the two bands is relatively strong. Stronger 
Coulomb repulsion in the tt band enhances the current 
density in the tt band, while it has little effect on the 
cr-band current density. 

Our model is suitable for describing MgB2. For pa- 
rameter values appropriate for MgB2 , we have found the 
above intriguing features in the LDOS and the current 
density, including additional bound states near the gap 
edge in the a band. The gap-edge bound states should be 
affected only weakly by the strong electron-phonon inter- 
action, as the energy of relevant phonons is much higher 
(above 60-70 meV) (Refs. [llllig^. Our predictions on 
the spectral properties of the a band can be tested by 
tunneling electrons onto the ab plane^i^ 
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